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The  purpose  of  this  paper  is  to  document  this  recent  emergence  of  fjeneralized 
tJiiC’ rkt'  as  a fundamental  computer-based  planning  tool  and  to  demonstrate  the 
power  of  the  associated  modeling  and  solution  technologies  when  used  in  aoncevt 
to  solve  real-world  applications.  To  accomplish  this,  the  paper  is  divided  into 
two  parts.  Part  1 focuses  on  how  generalized  networks  have  been  and  may  be  used 
to  model  a diverse  collection  of  significant  practical  problems.  First  we  discuss 
(in  nontechnical  terms)  the  model  structure  of  a generalized  network  (GN)  and  pro- 
vide a brief  historical  survey  of  applications  which  have  been  modeled  as  GN  prob- 
lems. Next  we  present  somewhat  newer  modeling  techniques  which  draw  heavily  on 
generalized  networks  as  a major  component.  The  pictorial  aspect  of  these  tech- 
niques has  proven  to  be  extremely  valuable  in  both  communicating  and  refining 
nonlinear  and  combinatorial  relationships. 

Part  II  discusses  the  design  and  analysis  of  large-scale  generalized  network 
computer  solution  techniques.  The  central  ideas  that  have  brought  GN  problems 
into  their  own  as  a fundamental  planning  tool,  by  providing  an  effective  method- 
ology for  computer  Implementations  are  clearly  outlined.  Part  II  further  contains 
an  in-depth  computational  study  of  solution  strategies  for  GN  problems  within  the 
framework  of  specializations  of  the  primal  simplex  method.  The  study  identifies 
an  efficient  solution  procedure  that  derives  from  an  integrated  system  of  start, 
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pivot,  and  degeneracy  rules.  The  resulting  method  is  shown  on  large  problems  to 

be  at  least  SO  times  faster  than  the  sophisticated  state-of-the-art  LP  system,  i 

APEX-III.  Finally,  the  stud>'  demonstrates  that  the  memory  requirements  of  the 

method,  as  well  as  its  solution  times,  are  sufficiently  small  to  warrant  its  use 

as  a computer- based  planning  tool  not  only  in  a batch  processing  environment,  but 

also  in  an  interactive  environment. 
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I nt  roiluct  ion 


.'hirijcfomont  Scionco  has  prop.rossfd  rapidly  since  World  War  II,  creating 
a virtual  explosion  of  new  knowledge  about  ways  to  solve  optimization  problems  in 
industry  ana  government.  The  single  mc'st  important  factor  motivating  this  explo- 
sion of  new  knowledge  has  been  the  parallel  growth  of  the  computer  industry.  The 
computer  has  given  management  scientists  the  capability  to  record  and  manipulate 
extremely  large  amounts  of  data  in  an  efficient  manner.  Without  this  capability, 
many  of  the  tools  of  management  science  would  be  mere  theoretical  niceties. 

The  advent  of  the  computer  has  given  rise  to  the  development  of  computer- 
based  planning  models.  The  techniques  for  building,  solving,  refining,  and 
analyzing  such  models  have  undergone  a steady  evolutionary  development  as  computer 
hardware  has  changed.  As  a result  of  this  evolution,  linear  programming  (LP)  has 
emerged  as  one  of  the  most  widely  used  tools  of  large-scale  computer-based  plan- 
ning. This  is  largely  due  to  the  fact  that  LP  models  portray  many  practical  ap- 
plications and  current  commercial  LP  computer  packages  have  been  able  to  solve 
problems  sufficiently  large  to  meet  the  requirements  of  real-world  settings. 

General  linear  programming  techniques  have  not  proven  to  be  totally  satis- 
factory, however.  The  conceptual  design,  data  collection,  and  computational 
capabilities  of  such  models  have  often  fallen  short  of  what  is  required  for  a 
truly  useful  decision  support  tool.  As  a result,  management  scientists  have  in- 
tensively studied  the  structures  of  LP  models  and  have  determined  fundamental 
building  blocks  which  appear  in  many  applications.  Network  models  have  long 
been  recognized  as  being  key  building  blocks  for  many  problem  formulations. 

Studies  of  the  structure  of  LP  models  from  practical  applications  have  underscored 
this  conclusion.  Because  of  this,  mathematical  programmers  have  devised  ingenious 
solution  algorithms  for  particular  types  of  network  problems.  The  computer  age 
has  spurred  the  development  of  computer  codes  embodying  these  algorithms  which 
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havi'  provfn  in  practice'  t<<  be  dramatically  more  ctfeoLivo  than  general  LI’  codes 

for  solving  these  problems  (1,  11,  19,  20,  28,  12,  18,  42  | . 

As  a consec|nenci' , various  subclasses  of  network  models,  such  as  shortest 
path,  maximum  flow,  minimum  spanning  tree,  assignment,  transportation,  and  trans- 
shipment models,  have  come  to  be  considered  as  rightfully  belonging  to  the  class 

Ilf  fundanu'ntal  computer-based  planning  tools.  That  is  to  say,  a model  is  cus- 

tomarily regarded  as  a fundamental  computer-based  planning  tool  if  it  is  capable 
of  accurately  depicting  a variety  of  significant  applications  either  directly,  or 
.IS  an  integral  component.  Additionally,  if  the  underlying  model  is  a special 
class  of  an  LP  model,  then  large-scale  algorithms  and  computer  codes  must  exist 
which  are  at  least  an  order  of  magnitude  faster  than  state-of-the-art  commercial 
LP  pack.'iges  for  this  model  type. 

Previously,  generalized  network  (CN)  models  [12,  14,  30]  have  not  been  con- 

sidered to  be  fundamental  computer-based  planning  tools.  Bradley  [10]  in  1975 
stated  that  GN  problems  ''in  the  near  future  . . . could  come  to  be  regarded  as  a 

r 

fundamental  model."  The  primary  reason  this  important  problem  class  has  not  pre- 
viously achieved  this  distinction  has  been  the  lack  v.f  efficient  large-scale 
computer  codes.  At  an  early  stage  in  the  evolution  of  the  computer  age,  research- 
ers developed  ingenious  special  purpose  solution  algorithms  for  these  problems. 
Further,  in  the  50's  the  early  pioneers  of  management  science  identified  important 
problem  classes  which  could  be  modeled  as  GN  problems.  Unfortunately,  these  model 
Identification  efforts  did  not  continue  with  any  rigor  in  the  60's  because  no 
effective  special  purpose  computer  codes  appeared.  Thus,  such  models  could  only 
be  solved  with  general  purpose  LP  codes. 

Recently,  however,  the  situation  has  been  changing  and  the  somewhat  overdue 
computer  codes  for  generalized  networks  are  at  last  appearing.  As  a consequence. 


miuiolers  have  bc'giin  to  devote  attention  to  delermininK  if  an  I,P  model  is  a GN 
problem  and,  more  importantly,  to  devising  formulations  in  which  generalized 
networks  play  the  role  of  critical  components. 

Tlie  purpose  of  this  paper  is  to  document  this  recent  emergence  of  generalized 
networks  as  a fundamental  computer-based  planning  tool  and  to  demonstrate  the 
power  of  the  associated  modeling  and  solution  technologies  when  used  in  concert 
to  solve  real-world  applications.  To  accomplish  this,  the  paper  is  divided  into 
two  parts.  Part  I focuses  on  how  generalized  networks  have  been  and  may  be  used 
to  model  a diverse  collection  of  significant  practical  problems.  First  we  discuss 
(in  nontechnical  terms)  the  model  structure  of  a generalized  network  and  provide  a 
brief  historical  survey  of  applications  which  have  been  modeled  as  GN  problems. 

Next  we  present  somewhat  newer  modeling  techniques  which  draw  heavily  on  generalized 
networks  as  a major  component.  This  collection  of  modeling  techniques  is  called 
the  NETFORM  (network  formulation)  concept  or  approach.  The  pictorial  aspect  of 
this  approach  has  proven  to  be  extremely  valuable  in  both  communicating  and  re- 
fining nonlinear  and  combinatorial  relationships.  Additionally,  the  NETFORM  concept 
often  yields  a formulation  that  enables  the  problem  to  be  solved  as  a sequence  of 
GN  problems  with  dramatic  gains  in  efficiency  over  alternative  approaches.  To 
provide  an  understanding  of  this  approach  and  the  role  of  generalized  networks 
within  it,  two  real-world  applications  are  described  which  have  profited  by  its 
use.  Thus,  Part  1 illustrates  the  manner  in  which  generalized  networks  are  emerging 
as  fundamental  building  blocks  for  modeling,  communicating  and  solving  a multitude 
of  problems. 

Part  II  discusses  the  design  and  analysis  of  large-scale  generalized  network 
computer  solution  techniques.  The  central  ideas  that  have  at  last  brought  GN 
problems  into  their  own  as  a fundamental  planning  tool,  by  providing  an  effective 
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mothodo lo^y  for  computer  imp lementat ions,  are  clearly  outlined.  i’art  II  further 
contains  an  in-depth  computational  study  of  solution  strategies  for  ON  problems 
within  the  framework  of  specializations  of  the  primal  simplex  method.  The  study 
identifies  an  efficient  solution  procedure  that  derives  from  an  integrated  system 
of  start,  pivot,  and  degeneracy  rules.  The  resulting  method  is  shown  on  large 
problems  to  be  at  least  SO  times  faster  tlian  the  sophisticated  state-of-the-art  LP 
system,  APKX-III.  (Thus,  the  metliod  can  solve  a problem  every  week  for  a year  and 
consume  the  same  amount  of  computer  time  required  to  solve  the  problem  only  once 
with  the  I.P  system.)  Finally,  the  study  demonstrates  that  the  memory  requirements 
of  the  metliod,  as  well  as  its  solution  times,  are  sufficiently  small  to  warrant 
its  use  as  a computer-based  planninq  tool  not  only  in  a batch  processing  environ- 
ment, but  also  in  an  interactive  envi ronment . 

PART  1 - GENERALI  7.F.D  NETWORK  MODEL 

1 . 0 PROBLEM  DEFINITION 

The  generalized  network  problem  represents  a large  class  of  LP  problems. 

This  class  Includes  any  LP  problem  whose  coefficient  matrix,  ignoring  simple 
upper  bound  constraints,  contains  at  most  two  non-zero  entries  in  each  column. 

A large  portion  of  the  literature  on  LP  problems  has  been  devoted  to  the  special 
cases  of  the  GN  problem  In  which  the  non-zero  elements  of  a column  consist  of  a 
1 and  a -1  (either  initially  or  by  linear  transformation).  This  condition  iden- 
tifies the  problem  as  a pure  network,  whose  Instances  Include  shortest  path, 
maximum  flow,  assignment,  transportation,  and  transshipment  problems.  The  GN  pro- 
blem, by  allowing  other  non-zero  doubletons  (and  singletons)  in  a column,  is 
actually  the  broadest  classification  of  linear  network  related  problems.  Practi- 
cal settings  In  which  such  GN  problems  arise  include  resource  allocation,  produc- 
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tiuii,  distribution,  sclit'dul  ins , capital  budRctins,  and  many  other  problem  types 
wbicli  will  he  elaborated  subsequently. 

The  most  effective  procedures  for  modeling  and  communicating  pure  network 
problems  are  based  on  viewing  these  problems  as  directed  graphs.  A generalized 
network  can  also  be  represented  as  a directed  graph  as  follows. 

I'nder  the  assumed  existence  of  a finite  optimum,  it  is  possible  to  transform 
the  coefficient  matrix  (by  scaling  or  by  complementing  a variable  relative  to 
its  upper  bound),  so  that  if  a column  has  two  non-zero  entries,  at  least  one  of 
these  is  -1.  In  this  way,  a directed  arc  is  "formed"  from  the  node  associated  with 
the  -1  to  the  node  associated  with  the  other  non-zero  entry.  (If  both  entries  are 
-1 , the  arc  may  be  directed  either  way.)  Columns  with  single  non-zero  entries 
give  rise  to  arcs  incident  on  only  one  node. 

There  is  an  important  distinction  between  arcs  in  pure  network  problems  and 
arcs  in  CN  problems.  An  arc  of  a generalized  network  has  a multiplier  associated 
with  it.  This  multiplier  is  the  non-zero  coefficient  associated  with  the  node  at 
the  head  of  the  arc  (i.e.,  to  which  the  arc  is  directed).  In  pure  networks,  this 
multiplier  Is  always  +1. 

Consider  the  following  GN  problem; 
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The  network  .issoc  iated  witli  this  problem  is  shown  in  Fisiire  1.  As  with  pure 
network  problems,  each  row  of  the  coefficient  matrix  is  associated  with  a node 
and  each  column  with  an  arc.  (In  other  words,  a node  corresponds  to  a problem 
equation  and  an  arc  corresponds  to  a problem  variable.)  The  arc  is  directed 
from  the  node  assoc  i.ated  with  the  -1  entry  toward  the  node  associated  with  the 
other  non-zero  entry.  Likewise,  each  arc  has  a cost,  lower  bound,  and  upper 
bound  which  are  shown  in  Figure  1 within  the  square  and  parentheses,  respectively. 
The  non-zero  multiplier  associated  with  an  arc  is  shown  in  Figure  1 within  a 
triangle.  The  constant  terms  (right  hand  sides)  of  the  problem  equations  iden- 
tify supply  and  demand  requirements  attached  to  the  corresponding  nodes.  A 
negative  constant  term  identifies  a supply  (which  by  convention  equals  the  ab- 
solute value  of  this  term),  a positive  constant  term  identifies  a demand,  and  a 
0 constant  term  identifies  a "conservation  condition"  in  which  the  amount  of  flow 
entering  the  node  must  be  exactly  matched  by  the  amount  of  flow  leaving  the  node. 

As  flow  passes  across  an  arc  in  a generalized  network  problem,  it  is  acted 
upon  by  the  non-zero  multiplier.  Thus,  the  amount  starting  out  on  an  arc  will 
not  necessarily  be  the  amount  arriving  at  the  opposite  end.  In  particular,  the 
multiplier  indicates  that  the  flow  entering  the  arc  is  multiplied  by  the  value  of 
the  multiplier  as  that  flow  leaves  the  arc.  For  example,  if  2 units  start  on  the 
arc  from  node  1 to  node  2 in  Figure  1,  4 units  will  arrive  at  node  2 since  the 
multiplier  is  2.  Likewise,  10  units  starting  on  the  arc  from  node  2 to  node  4 
will  result  In  2 units  arriving  at  node  4 since  the  multiplier  in  this  case  is  1/5. 
It  should  be  noted  that  the  arc's  cost,  lower  bound,  and  upper  bound  apply  only  to 
the  units  of  flow  enter iny  t/y*  arc. 
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2 . 0 APl’L  1 CAT  IONS  m OKNKKAl,  1 7^EI)  NKTWORKJi 

Oi'iu' r,i  I i /.cd  lU'tworks  can  be  used  to  model  tuimerous  |)roblems  for  wlilcb  there 
is  no  pure  network  etiuiv<ilent  . There  are  essentially  two  ways  in  which  the  mul- 
tipliers on  the  arcs  of  generalized  networks  may  be  interpreted.  First,  tney  can 
he  viewed  as  simply  modifying  the  amount  of  flow  of  particular  goods.  In  this  way, 
generalized  networks  can  represent  such  situations  as  evaporation,  seepage,  deteri- 
or.it  iviti,  breeding,  interest  rates,  sewage  treatment,  purification  processes  with 
varying  efficiencies,  machine  efficiencies  and  structural  strength  design.  How- 
ever, it  is  also  possible  to  interpret  the  multiplication  process  as  a transforma- 
tion from  one  type  of  good  to  another.  With  this  interpretation,  it  is  possible 
to  model  such  processes  as  manufacturing,  production,  fuel  to  energy  conversions, 
blending,  crew  scheduling,  manpower  to  joh  requirements,  and  currency  exchanges. 

The  following  applications  lend  insight  into  the  possible  uses  of  generalized  net- 
works . 

A complete  water  distribution  system  with  losses  has  been  modeled  by  Bhaumik 
[8]  as  a generalized  network  problem.  This  model  is  primarily  concerned  with  the 
movement  of  water  through  canals  to  various  reservoirs.  However,  the  model  must 
also  consider  the  retention  of  water  over  several  time  periods.  The  multipliers 
in  this  case  represent  the  loss  effect  due  to  both  evaporation  and  seepage. 

Turner  and  Gilliam  [17]  have  proposed  a file  reduction  model  which  has  the 
form  of  a generalized  transportation  model  with  a single  extra  constraint.  This 
model  is  designed  to  facilitate  the  reduction  of  extremely  large  microdata  files  to 
smaller,  statistically  representative  files.  The  objective,  in  this  case,  is  to 
minimize  the  amount  of  Information  lost  in  the  reduction  process.  The  arcs  repre- 
sent paths  from  the  original  records  to  the  reduced  records.  A non-zero  flow  on 
an  arc  implies  that  the  originating  record  Is  to  be  represented  by  the  terminal 
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reciird.  Tlic  multipliers  on  the  arcs  are  used  to  insure  that  the  reduced  file  is 
truly  representative  of  all  ot  the  original  records. 

Kim  [35 1 has  utilized  generalized  networks  to  represent  copper  refining  pro- 
cesses. The  electrolytic  refining  procedure,  in  this  case,  is  modeled  by  a large 
d-c  electrical  network.  The  arcs  are  current  paths  with  the  multipliers  repre- 
senting the  appropriate  resistances.  In  this  way,  Kim  analyzes  the  effect  of  short 
circuits  in  the  refining  process. 

Charnes  and  Cooper  |12)  identify  applications  of  generalized  networks  for 
both  plastic-limit  analysis  and  warehouse  funds-flow  models.  In  plastic-limit 
analysis,  tlie  network  is  generated  by  forming  the  equations  for  horizontal  and  ver- 
tical equilibrium  and  by  employing  a coupling  technique.  The  warehouse  funds-flow 
model  is  actually  a multi-time  period  model.  The  arcs  are  used  to  represent  sales, 
production,  and  the  inventory  holding  of  both  products  and  cash.  The  multipliers 
are  introduced  to  facilitate  the  conversions  between  cash  and  products. 

A cash  management  problem  has  been  modeled  as  a generalized  network  by 
Crum  [13].  This  model  for  the  multi-national  firm  incorporates  transfer  pricing, 
receivables  and  payables,  collections,  dividend  payments,  interest  payments,  royal- 
ties, and  management  fees.  The  arcs  represent  possible  cash  flow  patterns  and  the 
multipliers  represent  costs,  savings,  liquidity  changes,  and  exchange  rates.  Other 
applications  include  machine  loading  problems  [12,  lA,  A4],  blending  problems  [12, 
44],  the  caterer  problem  [14,  44],  and  scheduling  problems  such  as  production  and 
distribution  problems,  crew  scheduling,  aircraft  scheduling,  and  manpower  training 
[12,  14,  44]. 

3.0  INTEGER  GENERALIZED  NETWORKS 

The  uses  of  arc  multipliers  just  discussed  do  not  by  any  means  exhaust  their 
range  of  application.  In  fact,  upon  adding  the  requirement  of  discreteness,  which 
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lorcos  tin-  tliiws  (la  p.ifticular  arcs  Lo  occar  in  ialepcr  quantities,  the  flN' 
prohlem  is  capable  of  modeling  an  unexpected  diversity  of  problems.  For  example, 
iatroduciiiK  discreteness  into  the  (IN  model  produces  a framework  for  problems  such 
.is  scheduling;  v.iriable  length  television  commercials  into  time  slots,  assigning 
jobs  to  computers  in  computer  networks,  scheduling  pavments  on  accounts  where  con- 
tr.ictu.al  .agreements  specify  "lump  sum"  payments,  and  designing  communication  net- 
works with  capacity  c-onst  ra  ints . While  t Itese  are  "direct"  applications,  the  use 
ot  special  modeling  principles,  sometimes  c.ilied  NETFORM  principles,  enable  even 
somewh.it  more  complex  applications  to  be  modeled  and  solved  as  integer  GN  problems, 
in  fact,  the  NETFORM  principles  and  techniques  make  it  possible  to  model  any  0-1 
1.1'  problem  as  an  integer  GN  problem  (23,  27  ].  These  procedures  extend  quite  natur- 
al Iv  to  accommodate  mixed  integer  0-1  l.P  problems  where  the  continuous  part  of  the 
problem  is  a transportation,  transshipment  or  generalized  network  problem  itself. 
Reference  (A3J  illustrates  this  extension  and  shows  how  contemporary  financial 
c.ipital  allocation  problems  can  be  modelled  as  integer  GN  problems.  Many  other 
important  real-world  applications  have  such  a "mixed"  structure,  including  a 
variety  of  energy  models,  i>l;int  location  models,  and  phvsical  distribution  models. 

The  NETFORM  representation  is  mathem.it  ical  ly  equiv.ilent  to  any  of  the  alge- 
braic representations  that  c.in  be  arrived  at  by  customarv  mathematical  programming 
formulation  tecimiques . However,  by  the  creation  of  the  network-related  structures, 
which  may  or  may  not  be  implicit  in  any  customary  formulation,  the  NETFORM  repre- 
sentation permits  the  application  of  spec i .il i zed  solution  methods,  tailored  to 
employ  algorithmic  .advances  discussed  in  I’.irt  II  for  exploiting  the  graphical  re- 
lationships of  ni'twork  components.  The  remainder  of  this  section  briefly  describes 
the  b.isic  principles  of  the  NETFORM  appro, ich  ,ind  discusses  two  applications  which 


have  profited  from  its  use. 


Ki^Uirf  2 i 1 1 us  t r.'i  t fs  one  of  the  ustMul  modeling  devices  of  tlie  NKTFORM 
approach  that  finds  ,ipp 1 icat ion  in  a varielv  of  settings.  The  costs,  bounds, 
and  multipliers  are  represtuited  in  the  same  fashion  as  earlier.  In  addition, 
the  ustiTisk  on  the  .irc  from  node  0 to  node  A indicates  that  its  flow  must  be  an 
intei/(>r  amount.  Consequently,  in  view  of  the  upper  and  lower  bounds  on  this  arc, 
the  onlv  .acceptable  flow  v.ilues  are  ex.actlv  0 and  1,  and  the  multiplier  thus  en- 
sures th.it  either  0 or  T units  of  flow  are  transmitted  to  node  A.  Further,  the 
onlv  possible  way  to  distribute  3 units  of  flow  into  node  A is  to  send  exactly  one 
unit  to  e.ach  of  the  nodes  1,  2,  and  3 since  each  of  the  three  arcs  leaving  A has 
an  upper  bound  of  1.  Thus,  in  sum,  the  following  effect  has  been  achieved:  wh&n 

the  flow  on  the  arc  from  node  0 to  node  A is  0,  the  flow  on  each  of  the  three 
arcs  out  of  node  A is  0;  when  the  flow  on  the  arc  from  node  0 to  A is  1 , the  flow 
on  each  of  the  three  arcs  out  of  node  A is  1. 

It  should  be  noted  that  multipliers  m.ay  further  be  att.ached  to  the  arcs 
leaving  notie  A,  so  th.it  their  flows  may  be  further  transformed.  Thus,  for  example, 
the  flow  on  the  arc  from  node  0 to  node  A can  represent  an  investment  decision 
(invest  il  flow  = 1,  do  not  invest  if  flow  = 0),  and  the  flows  on  arcs  out  of  A 
c.in  represent  components  of  the  Investment  (e.g.,  particular  stocks  in  a portfolio, 
tr.icts  of  land  in  .i  real  estate  venture,  items  of  equipment  in  a manufacturing 
oper.it  ion,  etc.).  .Multipliers  on  the  letter  arcs  would  then  express  the  number  of 
items  of  eat  h tif  these  investment  components  th;it  are  obtained  by  selecting  the 
main  investment.  (For  example,  a pa.ticular  equipment  investment  may  be  composed 
of  six  tmichines  of  type  1,  eight  machines  of  type  2,  and  so  forth.)  The  combina- 
tion of  arc  multipliers  and  the  0-1  Integer  restriction  gives  rise  to  what  is  called 
an  tnteqvr  'fvneralized  ne-twork  or  a 0-1  cfetierali  zed  network.  This  NF.TFORM  modeling 
tool  has  a variety  of  Important  uses,  as  demonstrated  more  concretely  by  the  fol- 
lowing two  real-world  applications. 
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riRiirc  2 

Generalized  Network  witii  Inti'ger  Flow  Restrictions 


Ihc  Air  Force  requires  Undergraduate  Flight  Training  (I'FT)  graduates  to  take 
advanced  flight  training  before  their  first  operational  assignment.  ITi  graduates 
must  also  take  from  one  to  four  survival  training  courses.  Each  individual  has  a 
different  background  and  therefore  may  require  a different  mix  of  these  courses. 
Furttiermore,  both  the  flight  and  the  survival  training  courses  are  offered  only  at 
certain  times,  are  subject  to  enrollment  limits,  and  generally  have  prerequisites. 

A set  of  feasible  course  scliedules  is  identified  for  each  I'FT  graduate,  and  given 
a "cost  rating”  by  a computer  program  called  tiie  UFT  Pipieline  Scheduling  Model. 
Feasibility  and  cost  considerations  depend  on  addltic'nal  factors  such  as  attendance 
requirements  at  Combat  Crew  Training  courses,  various  modes  of  transportation,  the 
number  of  dead  days  in  the  pipeline,  the  opportunity  for  the  FT  graduates  to  take 
leave  as  desired,  etc. 

The  objective  is  to  select  a particular  course  schedule  for  each  UFT  graduate 
so  that  the  complete  set  of  schedules  selected  will  satisfy  all  class  enrollment 
limits  and  yield  the  smallest  total  cost.  To  solve  this  problem,  tlie  personnel 
manager  in  the  Training  Pipeline  Management  Division  munuallii  assigns  each  graduate 
to  one  of  his  feasible  schedules,  trying  to  assure  that  all  enrollment  limits  will 
be  satisfied.  Clearly,  this  Is  a difficult  and  time-consuming  task  to  do  bv  hand; 
further,  the  total  cost  of  these  manual  assignments  may  he  tar  from  optimal. 
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In  snarrli  of  ;i  hot  lor  approach,  the  Air  Force  developt'd  an  integer  program- 
ming; lormnlation  fc>r  this  problem.  However,  the  IF  formulation  turned  out  to  be 
almost  ti'tally  resistant  to  solution.  Consequently,  we  reformulated  this  integer 
programming  problem  as  a 0-1  CN  problem  shown  in  Figure  3. 

1 he  elements  of  this  diagram  may  be  explained  as  follows.  The  node  Mi 

represents  the  i-th  man  and  has  a supplv  of  exactly  1.  Each  man  node  is  connected 

by  arcs  to  its  set  of  man/schedule  nodes.  These  connecting  man/schedule  arcs  have 

a multiplier  a.,  equal  to  the  number  of  classes  in  the  schedule  and  a cost  c.. 

i]  IT 

equal  to  the  cost  of  assigning  man  i to  his  j-th  schedule.  The  asterisk  again  indi- 
cates that  flow  must  be  integer-valued. 

The  arcs  emanating  from  a man/schedule  node  in  Figure  3 lead  to  the  individual 
classes  making  up  the  schedule.  Each  of  these  arcs  has  an  upper  bound  of  one. 

Thus,  if  a particular  schedule  is  "selected,"  then  every  class  in  the  schedule  is 
also  automatically  selected.  The  objective  is  to  pick  a schedule  for  each  man  that 
will  minimize  the  value  of  the  assignments  on  the  overall  program,  subject  to  the 
upper  and  lower  attendance  limits  for  each  class,  expressed  as  bounds  on  the  arcs 
from  class  nodes  to  the  sink  nodes  of  . igure  3.  All  arc  costs,  except  for  those 
attached  to  the  man/schedule  arcs,  are  thus  equal  to  0. 

The  UFT  problem  typically  involves  120  men,  200  classes,  and  460  schedules, 
resulting  in  a 0-1  formulation  with  520  cc)nstraints  and  460  0-1  variables.  The  0-1 
GN  formulation  involves  the  same  number  of  0-1  variables,  and  introduces  an  addi- 
tional 2,200  continuous  variables  (arcs)  and  780  nodes.  Wliile  this  might  seem  to 
represent  a fair  Increase  in  size,  viewed  from  an  LP  problem  context,  it  actually 
represents  a relatively  small  GN  problem.  This  0-1  GN  problem  was  solved  using  a 
specialized  branch  and  bound  procedure  with  GN  subproblems.  The  optimal  solution 
was  often  found  and  verified  after  only  30  seconds  and  in  some  cases  only  required 
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a total  solution  t imo  of  10  sor<>tuls  f)ti  .a  t'lH  0000.  Thus,  t lio  problem  was  trans- 
formed from  one  that  had  been  extremely  ditfienlt  to  solve  as  an  integer  program 
ti'  one  that  was  solved  easily  as  a N'KTKORM. 

Re t tie  1 i n£  ^ic  1 ea  r Reaeto r s 

A mixed  integer  programming  problem  for  determining  the  minimum  cost  refueling 
schedule  for  nuclear  reactors,  modeled  bv  Kaxmerskv  ['54],  has  similarlv  benefited 
by  the  use  of  the  NKTFORM  concept.  Although  the  mixed  integer  programming  formu- 
lation bore  no  apparent  connection  to  networks,  we  discovered  a way  to  express  the 
problem  by  a convenient  NKTFORM  representation  after  interacting  closely  with  Dr. 
Kazmersky.  This  representation  was  not  only  equivalent  to  the  original  formulation 
but  also  succeeded  in  reducing  its  size.  Ve  will  forego  the  details  of  the  trans- 
formation of  the  original  problem  to  a 0-1  flN  problem  because  the  steps  are  some- 
what intricate  and  the  original  formulation  bv  itself  consumes  more  than  twenty 
pages  of  [34].  However,  we  were  able  to  exploit  the  0-1  GN  formulation  by  develop- 
ing a branch  and  bound  solution  procedure  which  emploved  the  GN  optimization  pro- 
cedures discussed  in  Part  II  of  this  paper.  Using  this  computer  system,  four  ver- 
sions of  this  problem  were  solved  using  data  supplied  by  the  TVA.  The  first  three 
versions  required  half  an  hour  to  two  hours  to  solve  on  an  IBM  370/168  using  MPSX. 
By  contrast,  these  versions  were  easily  solved  in  10  to  20  minutes  using  the  0-1 
GN  formulation  and  the  specialized  solution  approach.  The  fourth  version,  which 
Involved  173  constraints,  126  zero-one  variables,  and  511  continuous  variables,  was 
by  far  the  most  difficult.  Again,  using  MPSX  on  an  IBM  370/168,  the  original  mixed 
Integer  formulation  was  run  for  eight  hours  and  then  taken  off  the  machine  to  avoid 
further  computer  run  costs.  The  best  (minimum  cost)  solution  obtained  had  an  ob- 
jective function  value  of  $136,173,440.  With  a 30  minute  time  limit  imposed  on 
the  0-1  GN  solution  effort,  a solution  was  obtained  that  was  more  than  $10,000,000 


ihi-apiT,  witli  an  obji'clive  function  value  nl  $I  2b  , 1 74 , 727  . Consequently,  this 
applicaticui  shows  tluit  problems  too  complex  lo  lie  solved  n[)tima!ly  (within  prac- 
tical time  limits)  bv  standard  approaches  may  be  helped  suhs t ant  i al I y hy  t he 
NKTKOKM  app roac h . 

4 . 0 1 Mr  1.1  CAT  IONS  KOR  JTIK  OK  CN^  MODjil.S 

The  foregoing  development  provides  two  powerful  incentives  for  adopting  a 
GN  formulation  where  npiiropr iate . One  is  the  ability  to  conceptualize  a GN  for- 
mulation graphically  rather  than  algebraically.  This  pictorial  aspect  has  highly 
beneficial  effects  for  teaching  people  how  to  formulate  their  problems  and  for 
communicating  mathematical  models  to  nonscient  i f ic  users.  The  other  major  incen- 
tive is  the  ability  to  solve  GN  problems — .and  by  extension,  a variety  of  problems 
with  GN  components — with  a remarkable  degree  of  efficiency.  Part  II  establishes 
the  underlying  support  for  this  second  incentive.  Namely,  it  presents  an  in-depth 
computational  analysis  of  algorithmic  rules  and  computer  implementation  procedures 
which  are  capable  of  solving  large-scale  GN  problems  well  above  an  order  of  magni- 
tude faster  than  state  of  the  art  I.P  codes. 

PART  II  - DKSIGN  AND  ANALYSIS  OF  LARGE-SCALE 
GENERALIZED  NETWORK  COMPUTER  PROCEDURES 

1.0  OVERVi™ 

This  part  presents  an  extensive  computer  analysis  of  algorithmic  rules  for 
GN  problems.  Computational  studies  of  pure  network  solution  procedures  have  done 
much  to  advance  the  state-of-the-art.  Excellent  testing  has  been  performed  on 
transport.'it  ion  [18,  20,  28,  16,  19,  44]  and  transshipment  [4  , 5,  11  , 19,  26,  33, 

37  , 42)  computer  code.s.  These  studies  h.ive  provided  critical  Insights  into  the 
best  methods  for  solving  such  problems  as  well  as  providing  benchmark  data  for 
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tuture  solution  procoduros . 

Howovor,  then*  li.ivo  i>i>on  no  in-depih  studios  to  dnto  roncorning  the  mucli 
broader  ilass  of  GN  problems.  This  is  not  to  say  tliat  computer  codes  do  not 
exist  for  solving  such  problems.  Code  development  iias  been  reported  by  l!isomann 
115],  Maurras  (40],  Glover,  Klingman,  and  Stutz  [25],  Bhaumik  and  Jensen  [9], 

1-anglev  [18],  and  Balachandran  [2],  among  others.  Most  of  these  papers  report 
findings  for  only  certain  classes  of  GN  problems  and  all  t>f  them  are  limited  in 
the  scope  of  the  computational  analysis.  Thus,  an  important  body  of  empirical 
research  is  lacking  in  the  network  literature. 

The  code  NETG  reported  by  Glover,  Klingman,  and  Stutz  [25]  was  selected  to 
form  the  basis  for  the  computational  testing  of  this  studv.  This  code  is  an 
implementation  of  the  highly  efficient  extended  augmented  predecessor  index  (EAPI) 
procedure  (18,  24],  and  embodies  many  of  the  latest  advances  in  solution  methodology 
for  generalized  network  problems.  In  addition,  the  code  is  written  entirely  in 
FORTRAN,  has  been  thoroughly  debugged  and  tested  during  more  than  three  years  of 
developmental  work,  and  is  organized  to  allow  parameter  manipulations  and  sub- 
routine insertions  to  be  carried  out  conveniently.  These  features  are  extremely 
important  to  an  empirical  study  of  strategic  Implementation  possibilities. 

In  any  computer  implementation,  there  are  numerous  steps  that  can  be  performed 
in  alternative  ways.  Experience  from  previous  studies  of  pure  network  problems  has 
shown  that  the  determination  of  an  effective  set  of  decision  rules  to  handle  such 
alternatives  can  have  an  enormous  Impact  on  the  efficiency  of  the  Implemented 
solution  method.  Consequently,  a primary  objective  of  this  study  is  to  Investigate 
decision  rules  for  the  GN  problem  and  establish  their  relative  merits.  The  rules 
determined  to  be  best  have  been  Integrated  to  produce  a code  which  has  been  tested 
against  the  highly  efficient  linear  programming  system,  APEX-III.  The  results  of 
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tins  c'oiTi]i;ir  i son  iiulioati'  tli.it  t lu'  sju'ci.il  puriioso  (;N  co<io  is  at  least  50  times 
faster  th.in  Al'EX-lll  on  l.ir>’e  ON  problems. 

2.0  1’ ROIlj.r^  STATEMENT 

We  m.iy  formalize  the  i 1 1 nst  r.it  ions  of  I’art  I liy  defininp  a fonora  1 i zed  network 
problem  ;is  follows: 


Minimi ze 

T 

C X 

(1) 

subject  to: 

(Jx  = b 

(2) 

0 < X < u. 

(1) 

where  each  column  of  the  coefficient  matrix  0 contains  at  most  two  non-zero 
entries.  It  will  be  assumed  that  u is  finite  (e.^t-,  usinp  "regularization"  (12) 
if  nece-ssary).  The  problem  as  stated  m.jy  be  conceptual  ized  as  a graph  containing 
vertices  .and  non-directed  edges.  However,  we  will  further  assume  that  if  a column 
of  0 has  two  non-zero  entries,  then  at  least  one  of  these  is  negative  and  by 
appropriate  scaling  has  a value  of  -1.  This  allows  the  problem  to  be  interpreted 
in  a directed  graph  (digraph)  setting,  as  noted  in  Part  1. 

Each  column  of  G is  associated  with  a directed  edge  (arc)  of  the  digraph  and 
each  row  is  associated  with  a vertex  (node).  An  arc  runs  from  the  taj^l  nod£,  asso- 
ciated with  the  -1  coefficient,  to  the  ^ead  node , associated  with  the  other  non- 
zero entry.  In  the  case  of  a single  non-zero  entry  in  the  column  of  G,  the  asso- 
ciated arc  is  called  a self-loop  and  is  incident  on  a single  node.  The  flow  on  the 
arc  is  defined  to  be  the  value  assigned  to  the  corresponding  component  of  x (i.e., 
to  the  variable  whose  column  of  G gives  rise  to  the  arc).  Thus,  by  this  association 
between  arcs  and  variables,  each  arc  has  an  associated  cost  per  unit  flow  (the 
corresponding  component  of  c)  and  an  upper  bound  on  the  flow  (the  corresponding 
component  of  u) . The  coefficient  in  the  column  of  G associated  with  the  head  of 
the  arc  is  called  the  arc  multiplier. 
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Tlif  r i^tit -h.iiicl  side  v.iliu*  ( coininiruiit  ot  t> ) assoc  i at  f<l  with  a particular 
node  dett-rraiiies  whetlu'r  flow  will  be  inserted  or  remove(i  t rom  the  network  at 
that  node.  If  the  right-hand  side  value  is  negative,  flow  is  inserted  and  the 
node  is  called  a source  node.  Likewise,  if  the  right-hand  side  value  is  positive, 
the  node  is  called  a .sink  node  and  flow  is  rem<ived  at  that  point.  If  a particular 
node  h.is  both  arcs  heading  into  it  as  well  as  out  of  it,  then  it  is  called  a 
t ran_s.stupmeju  ru)de.  All  other  nodes  will  be  either  pure  sources  or  pure  sinks. 

Since  the  computer  code  NKTC,  on  which  this  stud)  is  based,  employs  a 
specialization  of  the  primal  simplex  method,  a brief  discussion  of  this  speciali- 
zation is  in  order.  We  begin  by  examining  the  basis  structure  for  this  special  i- 
zat ion . 

It  may  be  assumed  without  loss  of  generality  that  C has  rank  m,  where  m is  the 
number  of  rows  in  G.  iftherwise,  the  problem  is  equivalent  to  a set  of  disjoint 
minimum  lost  flow  networks  (see  [21]).  Any  basis  for  the  problem,  then,  will  be 
composed  of  m linearly  independent  column  vectors  selected  from  G.  Graphically,  a 
working  basis  of  the  simplex  method  is  a set  of  m arcs  incident  on  the  m nodes  of 
the  problem.  It  has  been  shown  [16,  24,  38,  40|  that  the  graph  of  a basis  will  be 
composed  of  ;i  set  of  disjoint  quasi-trees.  Kach  quas i-t ree  is  a simple  tree  to 
which  a single  arc  is  added.  The  additional  arc  forms  exactly  one  cycle , which  is 
a unique  series  of  arcs  leading  from  a node  back  to  that  node.  By  convention,  to 
allow  this  characterization  to  apply  to  the  case  in  which  the  additional  arc  is  a 
self-loop,  a self-loop  is  regarded  as  a cycle. 

The  KAPI  method,  on  which  NETC  is  based,  is  specifically  designed  to  store  the 
bases  of  generalized  network  problems  in  the  graphical  quasi-tree  form.  All  of  the 
prinviil  simplex  operations  are  carried  out  by  working  with  the  basis  graphs  which 
are  stored  using  linked  list  procedures.  The  Inherent  advantages  of  this  method 
will  be  fully  described  subsequently. 
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3 . 0 D/Vl’^  I N THK  KAIM  MK  I HOD 

A gi'nora  1 i zt'd  network  i)roblem  is  simply  n type  of  I.P  problem  and  can  thus 
1)0  solved  using  any  standard  1.1’  solution  technique.  Improvemi‘nts  in  inversion 
and  reinversiun  processtus,  data  compact i f icat ion , and  pivot  selection  strategics 
have  provided  dram.itic  increases  in  the  efficiency  of  primal  simplex  computer 
codes  in  recent  years.  In  many  cases,  special  structure  such  as  embodied  in  the 
generalized  network  problem  is  detected  by  these  codes.  This  information  is  then 
used  to  reduce  stor.age  requirements  and  to  simplify  operations.  However,  none  of 
the  current  I.P  systems  is  capable  of  fully  exploiting  the  structure  of  generalized 
network  problems. 

07ie  of  the  conspicuously  exploitable  features  of  generalized  network  problems 
is  the  si>arsity  of  the  coefficient  matrix,  and  current  LP  codes  are  of  course 
designed  to  take  advantage  of  sparsity  to  store  data  economically.  If  the  problem 
has  been  transformed  to  digraph  form,  however,  storage  may  be  further  reduced.  By 
the  use  of  simple  ordered  lists  to  capture  the  digraph  structure,  NETG  is  designed 
to  store  only  the  cost  coefficient,  the  non-zero  multiplier,  and  the  upper  bound 
lor  each  column  of  the  coefficient  matrix.  In  this  way,  problem  data  can  often  be 
resident  in  fast  main  memory  even  for  extremely  large  problems. 

As  noted  in  section  2,  bases  for  generalized  network  problems  have  a very 
special  structure.  With  possible  reordering  of  the  rows  and  columns,  the  basis 
matrix  forms  a block  diagonal  matrix.  Each  of  the  blocks  is  either  triangular  or 
near-triangular  and  can  be  represented  as  a quasi-tree.  Johnson  [31,  32]  originally 
proposed  a linked  list  procedure  for  storing  simple  trees  and  suggested  its  use  for 
the  more  complex  quasi-trees.  Effective  labeling  procedures  for  restructuring 
(updating)  quasi-trees  by  reference  to  such  lists  are  the  core  of  the  EAPI  method 
developed  by  Glover,  Klingman,  and  .Stutz  1 24  j . 


Tlif  original  iiroblom  daL.i  (compac l 1 y storoci)  and  tlif  h.isis  matrix  (stored  via 
linked  lists)  are  the  only  data  elements  that  need  to  be  kept  by  a specialized 
code  suih  as  NKTC.  The  main  factor  here  is  that  a basis  inverse  need  not  be 
maintained.  Inverses  gener.illy  require  considerable  amounts  of  storage  and  in- 
volve numerous  (error-producing,)  arithmetic  operations  to  utilize  and  maintain  them 
Instead,  the  specialized  labeling  rules  of  the  KAIM  method  operate  on  the  basis 
graph  in  a manner  that  obviates  the  use  of  a basis  inverse. 

The  KAPl  method  124]  orients  each  tpjasi-tree  so  that  the  cycle  is  concep- 
tually located  at  the  top,  with  attached  trees  hanging  downward.  Tliis  orientation 
is  called  a rooted  cycle.  The  linked  lists  and  labeling  procedures  provide  the 
means  of  traversing  the  tree  in  both  upward  and  downward  directions. 

Without  detailing  minutely  the  rule.s  and  processes  of  these  procedures,  it 
is  useful  to  sketch  their  main  functions  in  overview.  In  particular,  the  two 
types  of  quasi-tree  traversal  accommodated  by  the  procedures  are  used  to  carry 
out  basis  exchange  steps.  The  first,  upward  traversal,  is  associated  with  opera- 
tions normally  requiring  pre-multiplication  by  the  inverse,  such  as  determining 
the  representation  of  an  entering  vector  (arc).  This  oi>eration  Is  performed  by 
traversing  the  unique  paths  from  selected  vertices  up  to  their  associated  cycle(s). 
Simultaneously,  the  equivalent  triangular  system  of  equations  is  solved,  in  effect, 
by  back  substitution.  A directed  trace  of  tiie  cycle(s)  using  the  cycle  factor(s) 

1 22,  24  I completes  the  operation.  Tiie  process  ot  determining  the  representation 
of  a vector  requires  traversing  at  most  two  quasi-trees  and  only  generates  the 
non-zero  entries  of  this  representation. 

Downward  quasi-tree  traversal  is  analogous  to  post-multiplication  by  the  in- 
verse. This  operation  is  used  to  calculate  updated  dual  variable  values.  In  a net 
work  interpretation,  there  Is  a dual  variable  associated  with  each  of  the  nodes  in 


till'  I'roblfin. 
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Tliosi'  .iri'  often  referred  to  in  the  literature  as  node  potentials. 

At  each  iteration,  new  dual  values  associated  with  .a  subset  of  the  nodes  in  a 
sin>;le  quasi-tree  must  be  determined.  A single  va  1 ui'  assori.itod  with  a node  on  the 
cycle  mav  be  determined  using  tlu-  cycle  f.ictor  |22,  24  | . The  resulting  system  of 
equations  is  triangular,  and  may  be  solved  by  traversing  downward  through  the 
quasi-tree,  again  employing  back  substitution.  This  process  (uily  involves  those 
dual  values  that  change  during  a basis  exchange. 

Therefore,  the  specialized  procedures  reduce  both  data  storage  requirements 
and  arithmetic  operaticnis  required  to  carry  out  basis  exchange  steps.  However, 
there  are  several  alternative  methods  for  implementing  the  graph  processing  steps 
that  underly  these  procedures.  The  primary  purpose  of  this  study  is  to  determine 
which  of  these  implementation  methods  art*  the  most  efficient. 

4 . 0 K.\ 

The  computer  code  NKTT;  is  coded  entirely  in  standard  FORTRAN  IV.  Tile  use  of 
machine  dependent  operations  is  totally  .avoided  in  order  to  ease  the  transition 
to  various  computers.  However,  this  fact  also  lends  credence  to  the  results  ob- 
tained from  different  computers,  since  none  has  a particul.ir  coding  advantage.  The 
program  was  initially  coded,  debugged,  .tnd  tested,  using  the  Rl'N  compiler  on  a 
CDC  bbOO  computer  vith  ,a  maximum  main  memory  allocation  of  110,000  words.  The 
.ompiete  capacitated  algorithm  occupies  8N  + 4A  + 8500  words  of  central  memory, 
where  N is  the  number  of  nodes  .and  A is  the  number  of  arcs  in  the  specific  problem 
being  solved. 

Since  most  of  the  testing  performed  is  of  a compar.it  ive  nature,  it  was  de- 
sirable to  obtain  a set  of  problems  that  met  certain  specifications  and  that  could 
be  made  available  on  a repeated  basis.  For  this  reason,  a generalized  network 
problem  generator  was  developed.  This  code  Is  a logical  extension  of  the  NKTGKN 


1 '7 ) prohli'm  >',i'iu’r.iLi>r  for  puro  network  prol)lomK.  NFTCIION  \ises  a random  number 
generator  to  oreati'  problems  that  eonform  to  ci'rtain  strin  tural  characteristics 
supplied  by  tlu'  user.  The  usi'r  can  tumtrol  the  number  of  nodes,  sources,  sinks, 
transshipment  nodes,  and  arcs  as  well  .as  the  cost  range,  the  total  supply,  atid  the 
I'apacitv  ri  s t r ic  t i ons  . Tlu-  tu-w  i-iuie,  N’F.TtlKNC. , has  the  added  feature  that  the  user 
can  specity  a range  ot  valiu-s  1 rom  whicir  the  arc  multiplier  values  are  chosen. 

Using  NKTUKNd,  the  problem  sit  in  I.ible  I was  generated  for  the  purpose  of 
com])utat ional  testing.  The  information  shown  is  the  input  specification  for 
NKT('d.NG.  It  should  be  possible  to  recreate  this  data  set  exactly  on  virtually  any 
medium  to  large  scale  computer.  The  stated  problems  were  specifically  chosen  so 
that  the  elfects  of  problem  structure  on  solution  time  could  be  noted.  The  pro- 
I bJem.s  v.iry  in  .sizi-  from  JOO  nodvs  and  ) 100  arcs  u;>  to  1000  nodes  and  7000  arcs. 

The  200,  500,  400,  and  1 OflO  node  problems  consist  of  different  types  of  capacitated 

and  unc.ipac  i t ated  gener.ilized  t r.inspor  t a t ion  and  transshipment  problems. 

With  the  problem  set  complete,  the  process  of  selecting  decision  rules  began. 
The  types  of  decision  rules  tested  iii'  1 ude  start  i>rocedures,  pivot  selection  tech- 
niques, methods  for  exploiting  degeneracy,  tolerance  levi-ls,  I5ig-M  values,  and 
I pivot  tie-breaking  rules.  It  was  econom  i c.i  1 1 v inte.-isihle  to  perform  tests  on  all 

possible  combinations  ot  dec  ision  rules  or  even  to  test  each  rule  on  all  of  the 
generated  problems.  For  this  reast>n,  it  was  imperative  tt>  design  the  tests  to 
yield  useful  and  valid  results  while  minimizing  computat ion;il  time.  An  outline 
of  each  test  is  given  in  the  next  section. 
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Tlu'  tfst  in>j  ot  al  tirnat  i VC  decision  rules  was  pertormed  on  a CDC  6600  com- 
puti.r  lo.ati'd  at  Liii'  University  of  Texas  at  Austin.  In  each  of  the  comparative 
tests,  an  attempt  was  made  to  execute  t'ne  codes  involved  during  comparable  time 
periods.  However,  it  should  be  noted  that  minor  differences  between  two  solution 
times  should  be  statistically  ignored.  Tor  this  reason,  many  tests  were  run  on 
large  problems  so  that  timing  variations  become  less  significant.  Kach  code 
version  makes  use  of  a real-time  clock  routine  supplied  by  CDC.  This  routine  can 
be  employed  using  a K0RTRi\N  subroutine  call  and  is  generally  accurate  to  two  decimal 
p 1 aces  . 

Start  Procedurt-s 

Ihe  first  phase  of  testing  involves  a comparison  of  three  different  start  pro- 
cedures. All  of  the  starts  tested  are  based  on  techniques  that  have  proved  effec- 
tive for  pure  network  problems.  The  first  of  these  is  the  artificial  start  proce- 
dure. It  attaches  an  artificial  arc  (self-loop)  to  every  node  in  the  problem.  These 
artificial  arcs  are  then  assigned  extremely  large  (Big-M)  cost  coefficients.  The 
artificial  method  is  the  fastest  start  procedure  to  execute,  but  it  does  not  yield 
an  advanced  initial  basis. 

The  second  method  tested  is  the  sequential  source  minimum  (SSM)  procedure. 

This  method  sequentially  examines  each  node  in  the  problem.  If  the  node  has  asso- 
ciated supply,  flow  is  assigned  to  the  arc  having  the  least  cost  that  leads  to  a 
node  with  positive  demand,  or,  if  none  exists,  to  a node  with  zero  demand.  The  flow 
is  set  equal  to  the  minimum  of  the  supply,  the  upper  bound  on  the  arc,  or  the  de- 
mand (if  non-zero).  If  the  flow  on  an  arc  is  set  equal  to  the  supply  or  the  demand, 
the  associated  node  is  eliminated  from  further  consideration.  The  number  of  passes 
made  through  the  nodes  is  at  the  discretion  of  the  user.  If  the  process  is  termin- 


, Ill’ll  bi'fori'  siipfiiy  .ind  ili'ni.mcl  li.ivt'  lia-n  fxh.iiis  t t-d , tin-ii  .irlilicial  arcs  are 
.ippi'ndoii . For  t lio  purposi’s  ol  testing,  the  mimher  ol  passes  was  set  to  1,  2, 
b,  and  exliaustive. 

Tlie  last  start  method  tested  is  ttie  exliaustive  node  supply.  This  method  is 
similar  to  the  seijuential  source  minimum  in  the  way  it  assigns  flow  to  arcs. 
However,  the  procedure  continues  to  assign  flow  out  ol  a particular  node  until  the 
supply  at  that  node  is  exhausted  or  until  no  further  arcs  exist.  At  that  point, 
the  next  node  with  supply  is  considered.  Upon  completion,  remaining  supply  and 
demand  are  met  hy  appending  artificial  arcs. 

Each  of  the  start  methods  described  above  was  tested  using  two  distinct  pivot 
selection  criteria.  These  are  the  node  first  negative  and  the  node  most  negative 
criteria.  Both  methods  are  based  on  examining  the  non-basic  arcs  leading  out  of  a 
given  node.  The  node  first  negative  method  selects  the  first  encountered  pivot 
eligible  arc  for  the  basis  exchange.  The  node  most  negative  method,  on  the  other 
hand,  selects  the  best  pivot  eligible  arc  (in  terms  of  the  magnitude  of  the  updated 
cost  coelficient)  from  the  arcs  out  of  the  node.  All  other  code  parameters  were 
held  constant  in  all  of  the  tests.  The  results  of  testing  on  the  first  14  test 
problems  can  be  found  in  Tables  II  and  III.  Regardless  of  pivot  criteria,  the 
exhaustive  pass  SSM  procedure  is  predominately  the  best  start  method  in  terms  of 
resultant  total  solution  time.  It  provides  a reasonable  trade-off  between  the  time 
spent  selecting  an  initial  basis  and  the  time  recovered  from  a reduced  number  of 
pivots.  As  compared  to  the  artificial  start  procedure,  there  are  cases  in  which 
the  exhaustive  pass  .S.SM  method  reduces  total  pivots  by  as  much  as  61%  and  total 
solution  time  by  as  much  as  35%. 


f’ivDt  Selection  Criteria 


Tables  II  and  111  can  also  be  used  to  pi’rforni  an  initial  analysis  of  pivot 
select  i(.'ii  criteria.  It  should  be  noted  that  the  total  solution  times  in  Table  III 
strictly  dominate  those  in  T'.ible  II.  The  obvious  indication  is  that  it  is  worth- 
while to  spend  time  searching  tor  the  "best”  pivot  elij;ible  arc.  Selecting  the 
"best"  out  of  a single  node,  as  the  most  neg.it ive  method  does,  can  reduce  total 
solution  time  by  as  much  as  48%.  For  this  reason,  additional  testing  was  conducted 
to  try  to  find  the  best  pivot  selection  criteria. 

An  obvious  clioice  for  a pivot  selection  strategy  is  the  candidate  list  method. 

An  S-R  candidate  list  procedure  employs  an  array  of  length  R.  This  list  contains 
the  pointers  to  pivot  eligible  arcs  selected  by  using  the  node  most  negative  pro- 
cedure R successive  times.  After  each  of  S pivots,  the  best  arc  that  is  still  pivot 
eligible  in  the  list  is  selected  to  enter  the  basis.  If  tnere  are  no  eligible  arcs 
on  tile  list  or  if  the  list  has  been  used  S times,  the  list  is  refilled  by  calling 
Che  node  most  negative  procedure  R more  times. 

A number  of  variations  of  this  method  have  been  computationally  tested.  The 
first,  G1  , employs  a 1-1  candidate  list.  This  is  equivalent  to  the  node  most 
negative  procedure  and  is  included  as  a comiiarative  device.  C2  initially  uses  an 
8-16  list.  However,  if  the  list  cannot  be  filled,  the  size  of  the  list  is  set  to 
the  number  found,  and  S is  set  equal  to  1/2R.  tl3  is  similar  to  G2  except  that  if 
the  list  cannot  be  filled,  it  is  reduced  in  size  to  a 1-1  list.  The  fourth  version 
searches  through  a maximum  of  100  nodes  to  fill  an  8-16  list.  If  the  list  cannot 
be  filled,  R is  set  to  the  number  found  and  S » 1/2R.  G5  is  the  same  as  G4  except 
that  200  nodes  are  searched  before  the  list  size  is  altered. 

These  five  codes  were  tested  using  the  1000  node  problems  (15-18).  A 2-pass 
sequential  source  minimum  start  was  used  and  <all  other  parameters  were  held  constant. 


io 


riu  nsults  .m.-  shown  in  T.ihlo  IV.  It  is  ovidL'iit  that  (Id  and  (I'i  clt'arly  dominate 
till'  otlier  codes,  with  (12  having  a slight  edge  ovir  (I'j.  Note  that  by  employing  a 
(.Miididati’  1 ist  , tot.il  solnl  ion  t ime  was  reduced  liy  as  much  as  407  over  a simple 
nodi'  lU'gative  criteria.  This  test  indicates  that  candidate  lists  should  he  used, 
but  it  doi'S  not  show  the  best  initial  list  size.  This  determination  will  be  made 
subsequent ly . 

Kl(w  Update  I’rmtedures 

ihe  initial  version  ol  the  computer  code,  C2,  often  performs  wasted  operations 
in  the  process  of  updating  the  flows  on  basic  arcs.  First,  it  has  no  check  routines 
for  identifying  a degenerate  jiivot  during  the  calculation  of  a representation. 

Thus,  unnecessary  representation  components  are  found  and  flows  are  modified  by  a 
zero  amount.  Gb  was  designed  to  eliminate  this  problem.  it  identifies  a degener- 
ate pivot  and  totally  skips  the  flow  update  procedures  wtienever  possible.  Another 
version,  (17,  was  created  to  allow  for  storing  the  representation  vector.  Both  G2 
and  (ib  recalculate  the  representation  during  flow  updates.  G7  maintains  an  extra 
array  to  store  this  information. 

A comparison  of  G2 , G(),  and  G7  on  the  1000  node  problems  is  given  in  Table  V. 

All  three  versions  used  a 2-pass  sequential  source  minimum  start,  an  8-lb  candidate  1 

i 

I 

list,  and  all  other  parameters  were  held  constant.  It  Is  obviously  advantageous  to  j 

check  for  degenerate  pivots  as  Gb  does.  Although  the  total  number  of  pivots  is  not 
necessarily  reduced,  the  total  solution  time  can  be  reduced  by  as  much  as  25%.  Note 
that  the  number  of  pivots  differs  due  to  the  fact  that  tolerances  for  zero  values 
are  used.  Thus,  Gb  selects  the  first  "zero"  arc  to  leave  the  basis  while  G2  picks 
the  "minimum  zero"  arc  to  leave  the  basis.  However,  the  tests  on  G7  indicate  that 
It  is  not  advantageous  to  go  one  step  further  and  keep  a representation  array.  Ap- 
parently the  overhead  of  maintaining  this  extra  information  Is  not  offset  by  a 


reduction  in  other  calculations. 


(i.iiid  iclatf  I.isi,  1’aranu‘t  f rs 


All  import. lilt  cons  idt’ r.i  t i on  tii.it  was  not  rosolvt-ci  by  Ltio  pivot  criteria  .section 
is  till-  determiii.it  ion  ol  the  proper  candidate  list  p.irameters.  Tliis  involves  set- 
tiiij;  liotii  tile  B i 7.1’  ot  the  list  and  the  nunituT  of  times  the  list  is  to  be  used  be- 
fore retillin^;  it.  Candidate  list  sizes  of  1-1,  5-10,  H-16,  and  15-25  have  all 
been  tested  on  tile  1000  node  problems  iisinj;  the  Cb  code.  A five  pass  SSM  start 
w.is  employed,  wliile  par.imeti’rs  other  tlian  list  size  were  lield  constant.  A list 
size  of  5-10  can  be  seen  in  Table  VI  to  domina'e  all  other  strategies.  It  is  in- 
teresting to  note  that  a larger  list  size  does  not  guarantee  a reduction  in  the 
total  number  of  pivots. 

Pivot  To 1 e ranee 

Tolerance  Ievel.s  are  used  in  computer  programming  to  avoid  the  problems  caused 
by  round-off  error.  These  tolerances  give  r.anges  about  values  within  which  equality 
is  assumed.  Such  tolerances  must  be  used  in  the  testing  of  updated  cost  coeffi- 
cients. Negative  coefficients  Indicate  pivot  eligible  arcs;  liowever,  very  small 
negative  numbers  may  actually  be  equivalent  to  zero. 

In  order  to  find  the  best  tolerance  level,  v.alues  of  0.000001,  0.01,  0.5,  and 
1.0  were  tested.  Tliese  results  can  be  found  in  Table  VII.  Tiie  06  code  was  used, 
employing  a 5-10  candidate  list  and  an  exhaustive  pass  SSM  start  procedure.  It  is 
extremely  interesting  to  note  the  varying  effects  that  tliese  tolerance  levels  have 
upon  solution  times.  Yet  all  of  these  codes  arrived  at  solutions  with  tlie  same 
optimum  objective  function  value.  The  value  of  the  tolerance  dramatically  changes 
the  choices  of  pivot  eligilile  arcs.  Kxtremely  small  values  force  the  code  to  per- 
form a series  of  disadvantageous  pivots.  ()n  the  other  hand,  large  tolerances  allow 
the  code  to  initially  overlook  advantageous  pivot  selections.  The  best  strategy  is 


to  select  a moderate  tolerance  value  of  0.01. 


liij'-M  V.iluf 

Ihf  lin.il  pa  r.inu' t r r valua  lo  t)f  Lfstial  is  tlu'  Ki^-M  va  I ut* . None  of  the  codes 
devi  loped  iniplenieiits  ,i  I'liase  1-1’h.ise  II  procedure.  Thert'fore,  .an  exact  value  of 
15ig-M  must  he  si  leitiui.  All  of  the  test  problems  have  cost  ranj^es  of  1-100.  Kx- 
perience  with  pure  network  computer  codes  indicates  that  Big-M  should  be  set  .as 
small  tis  possible  (still  insuring  feasibility).  The  Big-M  values  were  tested  using 
the  lb)  code,  with  an  exhaustive  pass  SSM  start,  a 5-10  candidate  list,  and  a pivot 
tolerance  of  0.000001.  The  test  results  are  shown  in  Table  VIII.  Values  of 
10000,  2000,  1000,  500,  250,  150,  and  100  were  all  tested.  If  100  Is  eliminated 
because  of  resultant  in f eas i b i 1 i t ies  , a value  of  150  clearly  dominates  all  otners 
in  terms  ol  tot.al  solution  time.  Note  that  in  problem  16,  the  total  solution  time 
was  reduci-d  by  over  42/J  simply  by  changing  the  Big-M  value  from  10000  to  150. 

Breaking  Pivot  Ties 

The  l.ist  decision  rule  to  be  tested  was  the  one  for  resolving  ties  in  the  test 
for  a minimum  rati<i.  The  previously  described  codes  h.ive  simply  selected  the  most 
convenient  (often  the  first  encountered)  minimum  ratic’.  A new  code  version,  G8,  has 
been  developed  to  test  a heuristic  rule  for  breaking  pivot  ties.  The  rule  specifies 
tliat  the  ratio  with  the  largest  denominator  be  selected  from  .imong  the  set  of  tied 
ratios.  Note  that  t le  resultant  code  cannot  take  advantage  of  degeneracy,  and, 
therefore,  is  a modification  of  G1  and  was  tested  .against  Gl. 

The  ri'sults  of  testing  G8  on  .all  of  the  generated  problems  are  shown  in  Table 
IX.  In  the  majority  of  cases,  G8  reduced  the  ti>tal  number  iif  pivots  performed.  How- 
ever, this  reduction  w;is  not  enough  to  offset  the  increased  number  of  operations 
Involved  witii  the  pivot  selection. 
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1 58.70 

3999 

(4 

it. 93 

4 j4  1 

1 51.51 

3570 

13 

1 39  4-) 

io5e 

1 59.46 

650 

1 1 

77.00 

73i3 

6o..O 

6868 

17 

61  .te 

6^24 

51  . 

0072 

la 

L“'‘_ 

5955 

80.06 

6051 

* r’ms  •ts^ac’  are  In  CP»J  5«condi. 


2 . i sons 

II  is  liitticult  to  tolly  assess  the  elfiriency  of  a method  unless  that  method 
is  pLued  ill  comparison  witli  other  widi*iy  aicepted  methods.  For  this  reason,  NKTG, 
.ilong  with  the  jiri'viously  determined  decision  rule.s,  was  compared  with  a number  of 
well-known  linear  programming  computer  codes.  These  codes  liave  been  divided  into 
two  categories:  1)  specialized  network  codes,  and  2)  standard  linear  programming 

com[uiter  p.ickages. 

I he  tirst  comparison  was  against  several  specialized  progr.ims  tor  solving  pure 
network  problems.  Pure  network  problems  are  specializations  of  ..^.'alized  net- 
works in  which  a’l  of  the  arc  multipliers  are  1.  The  pure  network  codes  are  all 
equipped  to  take  advantage  of  the  simplifications  that  arise  from  the  unit  multi- 
pliers. NKTi;,  on  the  other  hand,  contains  the  cumbersome  mechanisms  for  handling 
non-unit  multipliers  and  quasi-tree  structures.  This  comparison,  therefore,  indi- 
cates the  relative  "cost"  of  these  additional  mechanisms. 

The  specialized  codes  used  in  the  comparison  were  ARC-II  [5],  FNKT- 1 [19], 

I/O  I'NFT- 1 ['i'ij,  Bennington  [7],  Boeing,  General  Motors,  .SHAKK  (41  1,  and  SllPERK 
[4J.  The  Boeing  and  General  Motors  codes  were  developed  by  the  respective  corpor- 
ations. All  of  these  programs  are  coded  in  standard  FORTRAN-IV.  The  tests  were 
performed  using  the  RUN  compilir  on  a GDC  6600  computer.  Each  of  these  programs, 
with  the  exception  of  1/0  PNET-1,  maintains  all  instructions  and  data  in  central 
memory.  I/O  PNE'1-1  uses  problem  data  from  external  disk  storage  and  is,  therefore, 
capable  of  solving  much  larger  jtroblem.s  th.in  the  other  codes. 

PNKT- 1 , from  which  I/O  PNKT-1  was  formed,  is  one  of  the  fastest  known  special 
purpose  primal  simplex  codes.  It  emuloys  a list  structuring  procedure  which  is 
more  advanced  than  the  one  used  by  NETG  (see  (26|).  ARC-1 1,  which  is  faster  than 
PNET-I,  uses  sophistlc.ited  primal  labeling  techniques  |5|.  SUPERK,  Boeing,  General 
Motors,  and  SHARE  are  all  variations  of  the  out -of-kl 1 ter  algorithm. 
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Mil'  prohli'in  si'l  for  t hi'  oompar  i son  consisLs  of  T)  piiro  lu'lwork  problt-ms 
piiu' ra  I I'll  l>v  NKI'dl-.N  | I7|.  I'lii'  input  spi-c  i t i la  t i ons  for  tlii'si'  problems  tan  be 
found  in  | 17 ) . Problems  l-l  are  100  x 100  transportation  prolilems;  problems 
t)-IO  ari'  I ')()  X 1)0  transportation  problems;  problems  11-15  are  200  x 200  assign- 
ment probliTis;  .uid  problems  lb-15  are  400-1500  node  capacitated  and  uncapae itated 
t ranssli  i pment  prob  1 ems  . 

1 tie  results  of  the  tests  are  shown  in  Table  X.  ft  can  be  seen  that  NKTO  is 
approximately  three  times  as  slow  as  ARC-1  I wtiich  is  by  far  the  fastest  of  tiie 
group.  This  result  gives  some  indication  of  the  additional  calculations  required 
by  a generalized  network  code  over  a pure  network  code.  An  interesting  result  is 
tliat  the  NKTC  times  are  roughly  comparable  with  the  SUPERK  times.  This  is  impor- 
tant because  SUPERK  is  one  of  the  fastest  out-of-kilter  implementations. 

The  second  phase  of  comparison  involved  NETG  and  a state  of  the  art  linear 
programming  code  called  APEX-111.  APEX-1  I I is  maintained  by  CDC  and  is  operational 
on  all  CDC  bOOO  series  and  CYBER-70  .scries  computers.  The  purpose  of  this  test  was 
to  indicate  the  advantages  that  specialized  procedures  have  over  standard  approaches. 
However,  it  should  be  noted  that  APEX-111  employs  highly  sophisticated  techniques. 

The  two  codes  were  tested  using  seven  problems  generated  by  NETCENG.  NETGENG 
is  capable  of  randomly  generating  generalized  network  problems  of  virtually  any  size. 
The  input  specifications  for  the  selected  problems  can  be  found  in  Table  XI.  They 
range  in  size  from  a 50  x 50  generalized  transportation  problem  to  a 1000  node 
generalized  transshipment  problem. 

To  perform  the  comparison  between  NETG  and  APKX-III,  a CDC  CYBER-74  computer 
was  used  and  NETG  was  compiled  using  the  CDC  FTN  compiler.  The  results  can  be 
found  in  Table  XII.  The  basis  of  comparison  for  these  tests  was  a quantity  called 
an  SBU.  Each  procedure  incurs  SBU's  based  on  the  amount  of  CPU  seconds  used,  I/O 
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Table  X 

Solution  Tiin«‘s  (in  secontis) 


Prob lems 

NHTG 

ARC- I 1 

I/O  PNET-I 

PNET- I 

BENN 

SUPERK 

CM 

SHARE 

Boeing 

1 

2.11 

. 78 

1 . 75 

1.17 

20.25 

5.68 

46.25 

17.76 

30.25 

2 

2.61 

.89 

1 .96 

1.35 

24.36 

6.47 

63.30 

21.34 

21.59 

3 

3.44 

1.01 

2.13 

1.74 

34.56 

6.87 

105.72 

26.16 

31.47 

4 

3.52 

.95 

2.35 

1.47 

31.45 

6.57 

70.  74 

25.13 

36.47 

s 

3.81 

1.25 

2.92 

1.73 

52.10 

6.77 

90.10 

30.97 

47.73 

b 

5.67 

2.11 

4.50 

3.06 

61.00 

11.05 

92.32 

46.40 

46.64 

7 

7.17 

2.23 

5.20 

3.57 

DNR 

12.86 

157.31 

65.92 

113.12 

8 

8.87 

2.99 

7.18 

4.20 

DNR 

13.69 

160. 71 

81.00 

175.10 

9 

8.46 

2.99 

6.76 

4.35 

DNR 

13.40 

158.01 

81.21 

186.99 

10 

8.72 

4.02 

7.37 

5.57 

DNR 

14.13 

197.82 

84.24 

184.75 

11 

7.31 

1.92 

3.71 

3.12 

17.44 

6.44 

35.67 

19.93 

30.39 

12 

8.32 

2.36 

7.75 

4.48 

20.  31 

6.47 

28.43 

21.17 

22.08 

13 

7.67 

3.13 

8.58 

4.91 

24.92 

7.25 

31.39 

25.81 

20.02 

14 

9.55 

2.96 

8.46 

5.56 

27.40 

6.95 

18.62 

24.95 

23.11 

15 

10.97 

3. 12 

8.60 

5.91 

DNR 

7.56 

23.48 

27.05 

21.08 

16 

3.71 

1.38 

2.60 

2.15 

11.77 

5.27 

60.27 

21.51 

15.05 

17 

7.49 

1.87 

3.65 

2.90 

20.10 

8.36 

96.66 

32/40 

64.64 

18 

4.44 

1.26 

2.65 

1.70 

11.31 

5.13 

61.54 

20.06 

18.31 

19 

5.98 

1.72 

3.61 

2.40 

20.62 

8.49 

DNR 

31.75 

61.07 

20 

5.12 

1.28 

2.47 

2.47 

10.38 

4.69 

DNR 

18. 1 1 

25.72 

21 

7.23 

1.83 

4.19 

2.46 

20.35 

7.96 

DNR 

32.60 

61.39 

22 

4.61 

1.26 

2.60 

2.01 

9.97 

4.60 

DNR 

17.91 

24.84 

23 

7.12 

1.67 

4.08 

2.74 

19.81 

7.91 

DNR 

32.66 

67.96 

24 

6.12 

1.52 

3.60 

2.91 

11.71 

5.59 

DNR 

25.27 

21.57 

25 

6.90 

1.83 

5.49 

3.96 

18.27 

8.37 

DNR 

33.19 

48.40 

26 

4.99 

1.08 

2.66 

4.15 

11.38 

5.51 

DNR 

25.05 

19.  34 

27 

6.88 

1.62 

4.03 

4.31 

16.37 

7.50 

DNR 

30.45 

41.98 

28 

15.11 

4.40 

10.34 

5.67 

DNR 

13.91 

DNR 

53.87 

83.98 

29 

18.72 

4.87 

12.51 

6.55 

DNR 

14.51 

DNR 

52.55 

117.83 

30 

17.77 

4.88 

15.32 

8.10 

DNR 

16.00 

DNR 

61.33 

152.21 

31 

18.26 

5.68 

13.45 

8.48 

DNR 

17.05 

DNR 

61.33 

135.73 

32 

29.77 

7.42 

20.  33 

13.59 

DNR 

22.88 

DNR 

78.63 

553.93 

33 

26.64 

7.82 

20.17 

17.65 

DNR 

25.89 

DNR 

101.92 

210.14 

34 

32.24 

8.21 

19.42 

14.86 

DNR 

25.42 

DNR 

92.25 

248.16 

35 

37.47 

8.81 

26.16 

17.13 

DNR 

29.96 

DNR 

DNR 

DNR 

NA  - 

Code  and 

data  would 

not  fit  in 

104,000 

words  of 

memory. 

DNR  - Did  not  run. 


Table  XII 


NHG  vs.  APCX-I  I I 


NETG  1 

APlX-I I 1 

Prob 1 em 

SBU's^ 

Cost^ 

SBU ' s 

Cost 

1 

7.51 

$1.35 

62.65 

S 11.28 

2 

7..  ■ 

SI  .31 

80.93 

S 14.57 

3 

9.70 

SI  .75 

94.72 

S 17.05 

4 

16.65 

S3. 00 

453.02 

S 81.54 

5 

14.74 

S2.65 

742.61 

SI33.67 

6 

22.55 

S4.06 

1044.34 

$187.98* 

7 

50.22 

S9.04 

1633.64 

$294. 06‘ 

^CYBER-74  System  Billing  Unit. 

^ComputeU  at  SO. 18  per  SBU. 

**Stopped  after  10,000  iterations. 

Objective  Function  Value  = 25,337,282. 
Optimal  Objective  Function  Value  = 3,354,927 

'^Stopped  after  10,000  iterations. 

Objective  Function  Value  = 1,340,958,349. 
Optimal  Objective  Function  Value  = 3,964,490 


o|H*r.it  ions  porlurmed,  and  central  memory  used.  In  this  way,  SIUl's  may  he  used  to 
compute  a total  cost  for  the  job.  Cost  fiypires  have  been  included,  based  on  the 
lowest  cnc  price  per  SBU,  $0.18. 

The  results  are  quite  remarkable,  especially  when  the  dollar  char>>es  are 
compared.  NF.TC  is  in  some  cases  more  thaii  bO  times  more  efficient  than  AI’EX-III. 
Troblems  b and  7 had  to  be  prematurely  terminated  on  APEX-Ill  (after  10,000 
iterations)  due  to  the  exhorbitant  processing  costs  involved. 
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The  parpose  of  this  paper  is  to  docament  this  rc'cent  emertf(*nce  of  generalized 
networks  as  a fandamental  compater-based  planning  tool  and  to  di'monstrate  the 
power  of  the  associated  modeling  and  solation  technologies  when  ased  in  concert  to 
solve  real-world  applications.  ^ To  accomplish  this,  the  pape'r  is  divided  into  two 
parts.  -I’art  1 focuses  on  how  generalized  networks  have  been  and  may  be  ased  to 
modi'l  a diverse  collection  of  significant  practical  problems.  First  we  discuss  (in 
nontechnical  terms)  the*  model  structure  of  a generalized  network  ((IN)  and  provide  a 
brief  historical  survey  of  applications  which  have  been  modeled  as  (IN  problems. 

Next  we  present  somewhat  newer  modeling  techniques  which  draw  heavily  on 
generalized  networks  as  a major  component.  The  pictorial  aspect  of  these  techniques 
has  proven  to  be  extremely  valuable  in  both  communicating  and  refining  nonlinear 
and  combinatorial  relationships. 

f’art  II  discusses  the  design  and  analysis  of  large-scale  generalized  network 
computer  solution  techniques.  The  central  ideas  that  have  brought  GN  problems 
into  their  own  as  a fundamental  planning  tool,  by  providing  an  effective  methodology 
for  computer  implementations  are  clearly  outlined.  Part  II  further  contains  an  in- 
depth  computational  study  of  solution  strategies  for  (VN  problems  within  the 
framework  of  specializations  of  the  primal  simplex  mrHhod.  The  study  identifies 
an  efficient  solution  procedure  that  derives  from  an  integrated  system  of  start, 
pivot,  and  degeneracy  rules.  The  resulting  method  is  shown  on  large  problems 
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13.  Abstract  (Continued) 


to  he  at  least  50  times  fastf'r  than  the  sophisticated  state-of-the-art  LP  sy.stem, 
APi;\-IIl.  Finally,  the  study  dcnnonstrates  that  the  memory  requirements  of  the 
method,  as  well  as  its  solution  timf'S,  are  sufficiently  small  to  warrant  its  use 
as  a computer-based  planning  tool  not  only  in  a batch  processing  environment, 
hut  also  in  an  interactive'  environment. 


